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Studying the nature and consequences of electron-phonon interaction in manganites is an area of 
intense ongoing research. Here, in an attempt to model charge and orbital ordering in manganites 
displaying C-type antiferromaganetism, we study cooperative Jahn-Teller effect in two-band one¬ 
dimensional chains in the regimes of both strong and weak electron-phonon couplings. These chains 
exhibit orbital ferromagnetism with only d ,^2 orbitals being occupied. At strong coupling and in 
the antiadiabatic regime, using a controlled analytic nonperturbative treatment that accounts for 
the quantum nature of the phonons, we derive the effective polaronic Hamiltonians for a single 
chain as well as for interacting identically-long chains. Due to cooperative effects, these effective 
Hamiltonians manifest a dominant next-nearest-neighbor hopping compared to the usual nearest- 
neighbor hopping and a signihcantly enhanced nearest-neighbor repulsion. For densities up to half 
Hlling, upon tuning electron-phonon coupling, interacting-chain [single-chain] Jahn-Teller systems 
undergo quantum phase transition from a charge disordered state to a conducting charge-density- 
wave state characterized by a wavevector k = (7r,7r) [fc = tt]. On the other hand, up to half 
Hlling, a weak coupling analysis reveals a transition from a disordered state to an insulating charge- 
density-wave state with a wavevector that depends linearly on the density; the ordering is analyzed 
within a Peierls instability framework involving the dynamic noninteracting susceptibility at nesting 
wavevector and phonon frequency. Our analysis provides an opportunity to identify the regime of 
electron-phonon coupling in manganites through experimentally determining the charge-ordering 
wavevector. 

PACS numbers: 71.38.-k, 71.45.Lr, 75.47.Lx, 71.38.Ht 


I. INTRODUCTION 

Understanding the wide variety of exotic long-range 
orders in transition metal oxides (such as manganites) 
and designing artificial structures (such as heterostruc¬ 
tures, quantum wires, and quantum dots) using these 
materials is of immense fundamental interest and also 
of huge technological importance. Perovskite mangan¬ 
ites Ri_a;Aa;Mn 03 (R = La, Pr, Nd, etc., A = Sr, Ca) 
are systems containing Mn^’*’ ions with one electron in 
doubly degenerate Cg orbitals; this results in cooperative 
Jahn-Teller (CJT) distortion (on adjacent sites) which 
lifts the orbital degeneracy and produces a cooperative 
occupation of orbitals. Consequently, a strong interplay 
ensues between charge, spin, and orbital degrees of free¬ 
dom. This interplay leads to closely competing energy 
states with a variety of spin textures (such as metal¬ 
lic ferromagnets and A-, C-, CE-, or G-type antiferro- 
magnets), charge orders [such as density-dependent and 
density-independent charge density waves (CDWs)], and 
orbital orders (such as ferromagnetic and antiferromag¬ 
netic orbital density waves )Ai^. 

There have been numerous studies of the many-polaron 
effects produced by quantum phonons in a one-band one- 
dimensional Holstein model^ which is a simple case of 
noncooperative electron-phonon interaction^A^. How¬ 
ever, controlled mathematical modeling of CJT quantum 
systems that goes beyond modeling localized carriers^ 
has remained elusive (at least to our knowledge). In fact, 
a controlled analytic treatment of cooperative-breathing¬ 
mode effects in a single-band, many-polaron chain has 


been developed only recently^. The situation is also 
complicated because of lack of conclusive identification 
of the strength of the electron-phonon interaction in ox¬ 
ides such as manganites. Evidence of strong local Jahn- 
Teller distortions pointing to strong electron-phonon cou¬ 
pling has been reported in manganites by direct tech¬ 
niques such as EXAFSi^, pulsed neutron diffractioni^, 
or through direct evidence of orbital ordering using res¬ 
onant X-ray scatteringii. On the other hand, weak cou¬ 
pling was inferred in the overdoped regime (i.e., x > 0.5) 
using transmission electron microscop} d^d^ , orientation- 
dependent transport measurements^^, terahertz spec¬ 
troscopy (to measure conductivity and permittivity ) ^^ ’^^ , 
and coherent synchrotron source (to measure optical 
conductivity)^. 

In the overdoped regime (x > 0.5) of manganites such 
as Lai_a;Caa;Mn 03 , at low temperatures, charge order¬ 
ing results along with orbital ordering and antiferromag¬ 
netism. The nature of the charge order has been a sub¬ 
ject of much debate. Earlier on, charge order was de¬ 
scribed in terms of a strong-coupling picture of stripes 
of localized charges where the mixed-valence manganese 
ions split into Mn^^ ions and Mn^^ ions2ii2^. However, 
based on the findings in Refs. [26l - [^ a phenomeno¬ 
logical picture of an extended CDW (instead of local¬ 
ized charges) was developed in Ref. The experi¬ 

ments indicating weak couplingi^^— for 0.5 < x < 
0.75 in Lai_a;Caa;Mn 03 , also seem to agree with this 
picture of weakly-modulated CDW characterized by a 
concentration-dependent wavevector q = {1 — x)a* with 
a* being the reciprocal lattice vector—. Furthermore, 
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in another manganite (i.e., Lao. 2 Sro, 8 Mn 03 ), charge or¬ 
dering wavevector q = 0.2a* has been reported^!. As 
regards magnetic ordering, in the overdoped regimes 
of Lai-xSrjjMnOa (for 0.65 < x < 0.95}^“— and 
Lai_xCaa;Mn 03 (for 0.75 < x < 0.87)^^— C-type anti- 
ferromaganetism has also been reported. 

In this paper, compared to the phenomenological the¬ 
ory of Ref. I 29 I we attempt at a microscopic theory in the 
overdoped regime. To reduce the daunting complexity, 
we assume C-type antiferromagnetism and invoke only 
the charge and orbital degrees of freedom to model the 
remaining possibilities of charge and orbital orderings. 
We investigate the regimes of weak coupling and strong 
coupling and show how the regime of coupling can be as¬ 
certained experimentally. Since manganites with C-type 
antiferromaganetic ordering can be effectively considered 
as chains, we study the effect of electron-phonon coupling 
in a single chain and in interacting chains. At strong 
coupling, we find that the charge ordering wavevector is 
always q = 0.5a* irrespective of the value of the doping 
x; whereas at weak coupling, the ordering wavevector has 
the density-dependent form q = (1 — x)a*. Thus, based 
on the experimental determination of the charge order¬ 
ing wavevectors at various doping values x in C-type an- 
tiferromagnets, we believe that our theory offers an op¬ 
portunity for identifying the regime of electron-phonon 
coupling. 

The paper is organized as follows. We derive a gen¬ 
eral Hamiltonian for CJT effect in three dimensions in 
Sec. II. In Sec. Ill, we derive the effective Hamiltonian 
for CJT interaction in one dimension and study its nature 
at various electron-phonon couplings. In the subsequent 
Secs. IV and V, we analyze the CJT effect in interacting 
chains for the case of strong coupling in the antiadiabatic 
regime and for the case of weak coupling in the adiabatic 
regime, respectively. Before closing, we discuss the con¬ 
nection between our system and other systems in Sec. VI 
and present our conclusions in Sec. VH. 


II. GENERAL CJT HAMILTONIAN 

The general Hamiltonian for the CJT system in man¬ 
ganites can be written as = Ht + H^p -I- Hi , where Ht 
is the hopping term, H^p the electron-phonon-interaction 
term, and Hi the lattice term. We start with an over¬ 
complete basis = 3a;^—r^,'0y = —r'^ — 

which satisfies the relation tpx + 'fpy + 'fpz = 0. The basis 
state '02 corresponds to the orbital depicted in FigH] 
The hopping term can be expressed in the above basis 
as: 


+ H.C.], (1) 


where ^ are creation operators at 

the site {i,j,k) for dy 2 , and d^^ orbitals, respec¬ 
tively. The labelling indices i, j, and k run along the x-, 
y-, and z-axes, respectively. The electron-phonon inter¬ 
action term can be written as: 

Hep — ^ j 

Tly2.j^ j kQy-i^j^k z;i,j,k]: (2) 

where g is the electron-phonon coupling, M is the 
mass of an oxygen ion, ojq is the frequency of optical 
phonons, and n^ 2 (y 2 ^^ 2 yi^ y i. = dl 2 (y 2 _^ 2 ).jj ^da: 2 (y 2 .z 2 );i,i,fc 
are the number operators. Furthermore, Qx-,i,j,k: 
Qy\i,j,k and Qz-,i,j,k are defined in terms of the dis¬ 
placements \Uxii^j^k ^ ^y^i,j,k ^ —l,fc; 

Uz-,i,j,k & Uz-,i,j,k-i\ of oxygen ions around (and in the 
direction of) the dx 2 ^ dy 2 , and dz 2 orbitals, respectively, 
as follows. Q x',i,j,k — '^x’,i,j,k ^x',i—\^j^k: Qy,i,j,k — 

'^y,i^j,k '^y;i,j — l,ki and Qz;i,j,k — '^z;i,j,k ^ 2 ;i,j,/c— 1 - Here, 
besides considering the displacement of the ions, we also 
consider their kinetic energy, thereby invoking quantum 
nature of the phonons. Then, the lattice Hamiltonian is 
given by 


Hi = 


M 

T 

K 




i,j,k 


T o ^ ^, [^x-.i,j.k '^y\ij,k (^) 


i,j,k 


where Ux-yj^k, Uyyj^k-, and Uz;i,j,k are the time deriva¬ 
tives of the oxygen-ion displacements Ux-ij,k, Uy-yj^k, and 
Uz-,i,j,k, respectively. 

The usual orthogonal basis states 02,2 2 and ipz^ are 

related to the over-complete basis states 0x, ijjy, and tpz 
as follows: 


ljjx^-y^ = 
022 = 02 - 





( 4 ) 


From Eq. o we get, 

0x = ^(\/30x2-y2 - 0 ^ 2 ), 

■02/ = - ^(V302;2_y2 -h0^2), 

02 = 02= • 

( 5 ) 

Next, using Eq. ®, we express the general Hamilto¬ 
nian in the orthogonal basis tpx^-y^ and 0^2 as follows: 
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dx‘^-y'^\i,j,k 
t 




^T 'j 

j + l,A;’ —y'^;i,j-\-l,k/ 




(Js f) (i:;«J+«■<=■> -'E(«'k.,M«.‘'i.-..«.«) (; o) J+<«) 


Hep = - -guJo\/2MuJo 


E /^t \ f Qx;i,j,k “t“ Qy;i,j,k “t" z;i,j,k \/^Qx;i,j,k ~t“ ^\/^Qy]i,j,k\ f 

\ -^^iQ.■,^,J,k + ^/^Qv■,^,3,k ^Qx■,i,J,k + iQy■,^,J,k Jv 


‘,i,j,k 

'x^-y'^;i,j,k 


( 7 ) 


and Hi is again given by Eq. ([3]). Here, it should be mentioned that an expression for in an alternate basis has 
been derived in Ref. [H; however, these authors consider classical phonons. 

In the subsequent sections, we shall study special cases of our above general CJT Hamiltonian, namely, the single 
chain and the interacting chains. 


III. SINGLE CHAIN CJT MODEL AT STRONG 
COUPLING 


We consider a one-dimensional Jahn-Teller chain with 
cooperative electron-phonon interaction along the z- 
direction and non-cooperative electron-phonon inter¬ 
action (of the Holstein-type^i^) along the x- and y- 
directions as shown in Fig. [T](6). The lattice term given 
by Eq. can be written for this case as follows: 


H 


CJT 


M 


r-2 1-2 1-2 1-2 

— -IT ^[^x-Q.k + '^x-,l,k + ^y;0.fc + '^y,l,k 


K 


k 

+Wy;l.fc + ^z;fe]- 


( 8 ) 


We define Q^-^k — T Qy^k — ^y; 0 ,fc 5 

Qx;k = '^x\l,k '^x; 0 ,k^ and Qy^k = '^y;l,k '^y, 0 ,k and 

incorporate these definitions in Eq. ([8]) to obtain 


H 


CJT 


_ ^ V 


-,{Qx-k + Qyk} + jlQty-k + Qxy-k} + "^t^k 


/b-2 


K 




k ^ 


^IQx-k + Qyk} + -jiQty-k + Qxy-k} + 


( 9 ) 


where Q^y.^. = Qx-k ± Qy,k- For the present single chain 
case, Eqs. and o reduce to the following equations: 

+ H-c.), ( 10 ) 


r 


and 


ffCJT 

ep 


gujQy/^A'IujQ 

= - ^ {(Wz;fc - Uz-k-l 
k 1 


) ^ Qxy,k J" 


^^2 . ^ ' h 




^ *5xy;fe(^l2.fcda:2_j^2.fc -I- H.C.^ 


(11) 


Next, we note that the center-of-mass displacement terms 
and Q'yi. as well as the center-of-mass momentum 
terms and Q'y.^ of Eq. m do not couple to the 
electrons [as can be seen from Eqs. (HOD and (HU]. Hence, 
for our single chain case, Eq. ([H]) simplifies to be 


H 


CJT 


= E 


k ^ 


E 

k 

E 


-Miilk -b -Kulj, 


1 A/ • , , 1 A , , ' 

_n+2 j_ r)+^ 

2 4 ^xy,k ^24 


k 


1 M ■ 


IK 


_r^n-2 _|_ 

2 4 ^xy,k ^24 


( 12 ) 


The general Hamiltonian for the present single chain CJT 
case can be expressed as follows by adding Eqs. m, 
(fTTl) . and (IT^ : 


jjCJT ^ jj 


CJT 


H 


CJT 


H 


CJT 


(13) 


Next, by using the following second-quantized represen¬ 
tation of the various displacement operators: 


A 


a;;fc + az;fe 


+ bk 


J 


2^ 


Q xy,k 


Ck 


Wo 


2^ 


Wo 
























FIG. 1. (Color online) Depiction of (a) one-dimensional single-band chain with cooperative breathing mode; (b) cooperative 
Jahn-Teller chain involving d ^2 and d^ 2 _y 2 orbitals; and (c) interacting cooperative Jahn-Teller chains with d ^2 and d^ 2 _y 2 
orbitals. For simplicity, only d ^2 orbitals are displayed in (b) and (c). 


in the above Hamiltonian of Eq. m, we obtain 
^CJT ^ + H.c.) 


3^0 


k 




+ + bk){nz2-k + ^n^2_y2.^) 

a/3 

2 ~b ^k^(,^^2.}^^x^—y^]k ~b H.C.) 
-ba;o^(a^a^;fe-b6[6fc-bc|,Cfe), (14) 


where n^ 2 .fe = and nx 2 _y 2 .j, = 

We will now adapt the well-known Lang-Firsov 
transformation^ for the above Hamiltonian so that we 
can perform perturbation in the polaronic (Lang-Firsov) 
frame of reference. The transformed Hamiltonian is given 
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by ifC'JT ^ exp{S)H^-^^ eM-S) where 

S = g 'y ^ [(Q'z:fc ^z-k')ij)‘z‘^\k ^2^;fe+l) 

k 

+ “ 6fe)(«22;fe + 3na:2-y2;fc)]- (15) 

Here, in our modified Lang-Firsov transformation, it 
should be noted that we have included only the density 
terms and ignored the orbital-flip terms {d}^- 2 .i^dx' 2 -y‘^-,k 
and its Hermitian conjugate) appearing in the interaction 
part of the above equation (fHl) . This choice is dictated 
by mathematical expediency to arrive at an analytic ex¬ 
pression. Then, the Lang-Firsov transformed Hamilto¬ 
nian is given by = Hq -|- Hi where 


Ho = ujQ ^(a^Uzjfc + blbk + c\ck) 
k 

-I- n^ 2 _y 2 .^) 
k 

^0 ^ ^ ^z‘^;k^x‘^ — y‘^;k ~t~ ‘^9 ^0 ^ ^ '^z^;k'^z^;k+l 




(16) 


where the term 2g‘^ujQj2k''^z^-,ki^z^-,k+i arises because of 
the cooperative nature of the interaction; furthermore, 
the attractive interaction term —^g^ujQ nz2-knx'^-y'^-,k 
will be negated by a much larger repulsive Coulombic 
term Cl ^z^;fe^a: 2 -y 2 ;fe because of which no site can 
have both the orbitals occupied simultaneously. The re¬ 
maining term of is given by Hi = H( -|- Hi^ with 

Hi = - 1} + H.C.], (17) 

k 


where 7± = exp[± 5 ( 2 a 2 ;fe-a^;fe_i-a^;fe+i)±f (6fc-5fc+i)] 
and 


29" ^(cl -b cfc) 


k 

n t 


d^2.kdx^-y^-,k 


g9(“I;fc_i-“I;fc+^'()g-9(a*;fc-l-a*ifc-l-&fc) _|_ JJ 


.(18) 


Now, to perform perturbation theory, we note that the 
eigenstates of Hq are given by |n, m) = \n)ei'S)\m)ph with 
10,0) being the ground state. We consider the strong¬ 
coupling case g^ >> 1 and the antiadiabatic regime 
t/ujo < 1 ; consequently, the coefficients of the pertur¬ 
bation terms H^ and H^^ in Eqs. and dn , re¬ 
spectively, satisfy the conditions te~~^ « uiq and 

^gujoe~i^^ << wo- Now, the second-order perturba¬ 
tion term [obtained using Schrieffer-Wolff transformation 
as mentioned in Eq. ( 6 ) of Ref. 13 is expressed as 


^ 


{^\phHi \ TTl)ph |0)p/i 

j^ph jpph 

-C^O — -^m 


In Eq. CH), the contribution of cross terms involving 
Hi and Hi^ is zero because the phonons do not match; 
hence, we get 


m 

+ E 


^Iph^l \'^}ph{'^\phH-^ 


- E 

±jn J-''> 


ph 


{0\phH(^\m)ph{m\phH[^\0)ph 

Ef - EDf 


( 20 ) 


We will first evaluate the term involving Hi^ in the above 
equation. After some algebra, we get the following ex¬ 
pression: 


(Ojphgi |m)p^(m|p/ig/^|0)p/i 


pph p: 
£/q — Hj, 


iph 


^ ^ ( |^7T,2:2;fc -f Hx^ — y^-^k 27T,2;2./g77-3,2_y2.^ . (21) 


We note that the coefficients of the terms 7722 .j., nx 2 -y 2 .k, 
and 7722 . 773 , 2 _y 2 .fc in the above equation are much smaller 
than the coefficients of the same terms in Eq. (UHl) ; con¬ 
sequently, we ignore the contribution from Eq. m in 
the expression for the effective Hamiltonian of the CJT 
chain. 

Next, after performing some tedious algebra (using 
considerations similar to those in Ref. 3 , the effective 
Hamiltonian can be obtained as: 


ttCJT _ 
^eff - 




where 


Wo \ 4 


X E [42;fc-l(l “ 27722.fc)d22;fc-Hl + H.C. 


g^wo-f—e [ 4 , 1 , 1 , 

LJo \ 4 4 


E '^z^ik'n-z^ik+i^ 


( 22 ) 


G3(2,2,- 


= F3 (2,2,E +F2(2,2) + 2E2 ( 2 ,^ 


m 


(19) 


+ 2Fi (2) + Fi 


1 

4 


(23) 
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and 


effective Hamiltoniani^: 


G. ( 4,l,l,v,v 


H. 


CBM 

eff 


- ( 4,1,1, i + 2F4 (4,1,1, + 2F4 (4,1, i I 

+ F4 (^1,1, i + F 3 (4,1,1) + AFs (^4,1, i 
+ ^3 (4, 1) + 2F, (^1,1, + 2F, (1, i 1 

+ 2F2 (4,1) + 2F2 (^4, + F2 (1,1) + 4F2 (^1, i 

+ F 2 + Fi (4) + 2Fi (1) + 2i^i Q] , (24) 


with 


=-te 39 " +H.C.) 

k 

®®'G2(2, 2) ^ [ 42 .fc_i(l - 2 n, 2 .fc) 42 .fc+i + H.c. 


-e 

ujo 


+2 


g^uJo + —e-^s G3(4,1,1) 

Wq 


^ ^ 77-22./j.74^2. 


fe+1- 


(27) 


fc 


On using the relationship of Eq. (ESI) , the effective Hamil¬ 
tonian for cooperative-breathing-mode chain [given by 
Eq. (I?7l) ] simplifies to be 

HCBM 


te 39 ^ 


Fn{o>l, . . . , CXn) — ^ ^ ^ ^ 


mi — 1 rrij 


(Q!ig^)’"E..(ang^)”^" 

TTJi!... ninlimi -I- ... -b m„)' 


dz'^-k + H.c.) 

k 

9 

te® 


4g2a;o 


X] [ 42 ;fc-i(l - 2?422;fc)42.fc+i -b H.c. 


A general term of the form G„(ai, a 2 ,..., an) can be 
expressed as G„(ai, 02 ,..., an) = Fn{ai,a 2 , • ■ •, a„) -b 
YlkZiJ 2 cFk{aci,ac 2 ,---,ack) where the summation 
over c represents summing over all possible "Gm com¬ 
binations of m arguments chosen from the total set of n 
arguments {ai, a 2 ,..., an}- We then obtain the follow¬ 
ing useful relationship (derived in Appendix lAT) : 

Gn(ai, a 2 ,..., a„) 

n—l 

Fni^ai , a 2 ; ■ • ■ 7 an) F ^ ^ ^ ^ Fj^ (c^ci ? ^C 2 ? ■ * ■: ^cfc) 


2g^wo , t 


t Sg^wo 


=39- 


y^n^2.j,n^2. 


fe-i-i- 


(28) 


k—1 c 


for g 


’>>1 e^"=i“-9" 


= E 


d 


E n 2 ■ 

i=i 


(25) 


Then, on using the above approximation for 
Gn(ai, a 2 ,..., a„) at large g^, Eq. (|2^ simplifies 
as follows: 

HO/T 


J. -i^o2 

te 45 


4" H-C.) 
k 

4 te® 


17 g^wo 


^ ^ j^d^2.^_4(l 277-22;/„)(i22;7,;_|_l -b H.C. 

k 


2g‘^uJo 4 t 


(26) 


t ISg^wo 


In contrast to the effective Hamiltonian of Eq. 
the relatively simpler case of single-band cooperative 
breathing mode [shown in Fig. (Ua)] yields the following 


Comparing Eq. (l26ll with Eq. (l28)) , we find that the co¬ 
efficients of next-nearest-neighbor (NNN) hopping terms 
are approximately equal and so are the coefficients of the 
nearest-neighbor (NN) repulsion terms. Hence, both the 
CBM chain and the CJT chain should exhibit similar 
behavior up to half-filling. We diagonalize these effec¬ 
tive Hamiltonians using a modified Lanczos algorithm^ 
[with antiperiodic (periodic) boundary conditions for 
even (odd) number of fermions] and calculate the struc¬ 
ture factor S{k) at the ordering wavevector tt. Upon tun¬ 
ing the electron-phonon coupling g in the antiadiabatic 
regime, earlier we found that the CBM model undergoes 
a second-order quantum phase transition from a Lut- 
tinger liquid to a CDW state at strong couplingii. From 
Fig. [21 for both the CBM model and the CJT model, we 
see that the S{Tr) curves coincide more or less. Hence, up 
to half hlling, both the models exhibit similar CDW tran¬ 
sition upon tuning electron-phonon coupling g. However, 
above half filling in the CJT chain, particles will occupy 
the d^ 2 _y 2 orbitals of the remaining sub-lattice because 
there is no repulsion between electrons in the and 
dx 2 -y 2 orbitals on NN sites [see Eq. dUD]. Thus, in con¬ 
trast to the CBM model, the particle-hole symmetry is 
broken for the CJT model! 

In the subsequent sections, we shall study interacting 
CJT chains at couplings that are strong and weak. 


IV. INTERACTING IDENTICALLY-LONG CJT 
CHAINS AT STRONG COUPLING 

In the interacting CJT chains depicted in Fig. (He), 
we have C-type antiferromagnetism (i.e., each ferromag- 
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and 


j,k 


“t“ '^{Qx;j,k “t“ '^yj,k ;j; —1,/c) ^ d^2 


.(32) 


Next, we observe that the displacement operator 
does not couple to the electrons in the above Eqs. (1311) 
and (|32l) : therefore, we drop terms involving Q'^.j j. from 
Eq. (1501) . Then, the total Hamiltonian for the interaction 
CJT chains is given by 

^ICJT ^ jjICJT ^ jjICJT ^ ^ICJT ( 33 ) 


FIG. 2. (Color online) Structure factor SOk) showing CDW • i i i- i 

transition in CBM chain and CJT chain when ^ = 0.1 and representing the relevant displacement operators m 

for (a)iV = 16, = 4; and {b)N = 18, = 6^” ^^cond quantized form as 


netic chain is antiferromagnetically coupled to the adja¬ 
cent ones); consequently, there is no inter-chain hopping. 
In these chains, hopping of electrons occurs only between 
the ^22 orbitals along the z-axis. Furthermore, the chains 
are assumed to have cooperative electron-phonon cou¬ 
pling along the y- and the z-axes, whereas, along the 
x-axis the electron-phonon coupling is non-cooperative 
and is of the Holstein-type^i^. Up to half filling, only 
the d ^2 orbitals would be occupied so that the system 
lowers its energy through the mobile electrons. Here, we 
consider the interacting chains for only up to half filling; 
therefore, we neglect the occupation of orbitals in 

our analysis. Hence, the lattice term of Eq. m for the 
present case reduces to 


H, 


ICJT 

I 


M 


= tE[ 


+ + '^y,j,k + Wzp./c]-(29) 

j,k 


+ k + 


We define Q^-j ^ — Ux'i^j^k~^^x',o,j,k and Qx-j^k — ^x',i,j,k 
Ux-,o,j,k and recast Eq. (I55|) as follows: 


H, 


ICJT 

M 

" T 


Y 2 


^y;j,k ^z-,j,k 


j,k ^ 


.(30) 


For the present system of interacting identically-long 
CJT chains, Eqs. m and 0 simplify to be 

^ICJT = -tY (dl^:yk+ldz^-,3,k + H.C.) , (31) 

3:^ 


'^z;j,k — 


Qx;j,k — 


^z-,j,k + 


y/2A'IujQ 

^x-,j,k + ^x-,j,k 

\/2^Wo 


, '^y;j,k — 


^y;j,k + h-J,k 

y/2A'I(jjQ 


(34) 


and substituting in the total Hamiltonian of Eq. (I33L 
we obtain 

= -tYidY,k+idz^;j,k + H.c.) 


j,k 


- gt^oY 


^^z;j,k ^z;j,k)(j^z'^;j,k ^z'^;j,k-\-l) 


A^y.j.k + h-J,k)in^2.j f, - 
+ Wo Y^Aj,k^z:-j,k + bl.j f.byj^k + cl;yfcCx;j,fc), 

3,k 

(35) 

where n^ 2 .j i^ = ^.^ 22 .^ j,. To obtain an effec¬ 

tive Hamiltonian by performing perturbation theory, we 
device a relevant Lang-Firsov transformation for the 
above Hamiltonian. In the polaronic frame of refer¬ 
ence, _ exp(S')i7^‘"‘^^ exp(—S') where the anti- 

Hermitian operator S is identified as 


^ = -9Y1 


i^z;j,k ^z',j,k^{j^z'^-,j,k '^z'^;j,k-\-l^ 


3,k 


+ ^i^y;j,k byj^k)(,'n,z'^-,j,k kl^2.j^ii^) 

+ ^^A,j,k-Cx;yk)n.^;yk ■ 


( 36 ) 





























Subsequently, the Lang-Firsov transformed Hamiltonian 
is given by H^CJT ^ ffiCJT 


H, 


ICJT 

0 


— Wo + by-J^k^V,j,k + cl..j i.Cx-,j,k) 


j,k 


j,k 

, 13-2 • 


-te 4 3 ^( 42 .^-fc+idz 2 ;yfc+H.C.), 

j,k 


where terms ‘2g^^Qj2ik''^z^jkn'z^-jk+i and 

k‘^z^\j,k'>^z^■,]+!,k arise because of coopera¬ 
tive nature of the Jahn-Teller interaction along the 2 ;- 
and the y-directions, respectively. Furthermore, the 


perturbation is given by 




where 


= —te 4 ■ 


X 

j,k 

-bH.c. 


^z^-,j,k+l 


dz^.,,,k{rrr^" - 1 } 


(38) 


y 2 

^9 Wo / ^ fc 

j,k 

(37) 


T -j.k _ 

^ = exp 


zb y( 2 a 2 jj 0>z\j,k—l ^z;j,k+l) 


^ ^ (^3;j-l,fe+l "b byj-l^k byj^k+l) 


zb 


2V2 


(^^x;j,k ^x;j^k-{-l) 


(39) 


Since we consider strong coupling and antiadiabaticity, 
we get the condition te~~^ << wq. Then, similar to 
the single chain CJT case [on identifying the eigenstates 
of as |n, m) = \n)ei (S' |m)p/i with |0,0) being 

the ground state], the second-order perturbation term 
is given as: 

o-(2) {0\phH(^-^^\m)ph{m\phH(^J^\0)ph 

^ Jj^ph jpph ■ 


pph jp: 


After tedious algebra, we get the the following expression 
for the effective Hamiltonian for interacting CJT chains: 


J 


H. 


ICJT 

eff 


— te 4® fc -|-H.c.) 


j,k 


t -IS. 2 
- P 2 y 

UJo 


G5 ^2,2, i, ^ [42;j.fe-i(l “ 2 n^ 2 .^_fc)(l - n^2y_i_fc)(l - n;,2y+i_fc)42y fc+i -h H.c. 


+2 


y^WoH-e =3 Gg 4,1,1 


Wo 


i,k 

111111 
8’8’16’16’16’16 


^j.n22 .j--I- -y^wo ^ n22y-^n22y_|_i j,, (41) 
j,k j,k 


where G„(q;i, 0 : 2 , • ■ ■, Q:„) are the same as those defined in the previous section. For >> 1 (based on derivations in 
the Appendix [AT) . the above equation reduces to 


jjICJT 

^_ia(,2 = ~ + H.C.) 

j,k 


te-T34 


4 tes 
17 g'^ojo 


[<^l 2 .j^_i(l - 2n^2y_fc)(l - n^ 2 .j_i ,^){l - n^ 2 .j+i_fc) 42 y ^+1 -|- H.c. 

j,k 


2g‘^uJo 4 t 


t 13y^wo_ 


ISq'^ 1 y^Wo lSo2 \ 

e43 y n^2.yf.nz2.y^i + o~^e / Xz^q.fc^z^q-n.fc; 

j,k j,k 


(42) 


which is exactly the same as the single chain CJT re¬ 
sult [see Eq. (l26l) ] except for the inter-chain repulsion 
due to cooperative effects along the y-axis. Further¬ 
more, in the above equation, the occupancy-projection 


factors (1 — n^ 2 y-j-i fc) [appearing in the NNN hopping 
term] arise due to the fact that the reduced hopping in¬ 
tegral te~~^ is significantly smaller than the inter-chain 
2 

repulsion strength 
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FIG. 3. (Color online) Plots of (a) density-density inter-chain 
correlation function IF(l,m); (b) density-density intra-chain 
correlation function VF(0,m); (c) structure factor S{Tr,kz)', 
and (d) structure factor ^(O, kz) evaluated at adiabaticity 
t/iOQ = 0.1, system size = 12 and particle number Np = 4. 


Although in the present work we focus on interacting 
chains in two-dimensions only, effective Hamiltonian for 
interacting CJT chains in three-dimensions can be also 
obtained similarly. 

Using the modified Lanczos algorithm of Ref. HO, in 
order to study the nature of the system as the electron- 
phonon coupling g is tuned, we diagonalize the Hamilto¬ 
nian in Eq. ra and calculate quantities such as density- 
density correlation function, structure factor, fidelity, 
and fidelity susceptibility in the following subsections. 
Furthermore, a two-dimensional system with interacting 
CJT chains is taken to be mimicked by a two-chain CJT 
system with periodic boundary conditions. Although 
we perform our calculations for a conservative value of 
t/u!o = 0.1, our results are valid for fairly larger values of 
t/ujQ < 1 as demonstrated in the Appendix IbI 


A. Density-density correlation function, structure 
factor, and order parameter 

The two-point correlation function, for density fluctu¬ 
ations of electrons that are apart by a vector {l,m), is 
defined as 

for a filling fraction v = {n^ 2 .j ^.) = where N is the 
total number of sites in the system and Np the total num¬ 
ber of the particles in the system. The Fourier transform 



0 


2 3 4 0 1,2 3 

m k, 


FIG. 4. (Golor online) Plots of (a) inter-chain correla¬ 
tion function lU(l,m); (b) intra-chain correlation function 
tU(0,m); (c) structure factor S(Ti,kz)', and (d) structure fac¬ 
tor 5(0, kz) at t/cjo = 0.1, A = 16 and Np = 4. 


0.2 

0.1 

I 


- 0.1 



FIG. 5. (Color online) Plots of (a) correlation function 
lU(l,m); (b) correlation function W{Q,m)\ (c) structure fac¬ 
tor S{'K,kz); and (d) structure factor S{0,kz) at t/uio = 0.1, 
A = 20 and Np = 4. 


of W(l,m), i.e., the structure factor S(ky,kz), is given 
by 

S{ky, kz) = m), (44) 

l,m 

where for our two-chain CJT system^ on assuming pe¬ 
riodic boundary conditions, ky = with 

Uy = 1, 2; and kz = with Uz = 1,2,y- 

We consider only even values of A/2 so as to obtain sub- 


A 

3 


■ w / u / \ r 

A 2 

' ^ 


' , ^ 

1 

n . 

Y/^ ^(c) 
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FIG. 6. (Color online) Display of (a) two-point correla¬ 
tion function (b) two-point correlation function 

1F(0, m); (c) structure factor S{TT,kz)', and (d) structure fac¬ 
tor S(0, kz) at t/uJo ~ 0.1, N = 24: and Np = 4. 




FIG. 8. (Color online) Order parameter S* (tt, tt) as a function 
of the coupling g when t/wo = 0.1. 


To begin with, we observe that 

S{tt,tt) = I X! “ ^ + W{0,m)]. 

V'^even ^odd/ 

Based on the derivation of Eq. (|C11I) in the Appendix [Cl 
we can write 

5(7^: - ^l.o) + (^2.0 - ^2.e)]"), (46) 

where e = Efc_„ {^ 3 ,o = is the 

number operator for the total number of particles in the 
even (odd) sites of the j-th chain. Next, again from Eq. 
(|45]) . we find 

SiO,TT)= ( ^ ^ ] [IE(l,m) + IE(2,m)]. 

\?7leven ^odd/ 

As shown by the derivation of Eq. (|C12ll in the Appendix 
[Cl 5'(0,7r) can be expressed as 

^(0,^) = - ^l.o) + {N2,e - N^,o)?). (47) 

Eurthermore, Eq. (HSl) also yields 


FIG. 7. (Color online) Depiction of (a) correlation function 
IT(l,m); (b) correlation function VF(0, m); (c) structure fac¬ 
tor S{n,kz)', and (d) structure factor S{0,kz) at t/u)o = 0.1, 
A = 24 and Np = 6. 


S{7T,0)= ^ ^ [-IE(l,m) + IE(2,m)]. 

\?71even ^odd/ 

Then, from Eq. (IC13|) in the Appendix [Cl we observe 
that 


lattice structure. We will now derive exact relations for 
the structure factor in a few special cases. To this end, 
the above equation is simplified to be 



[(-!)”“ W(l,m) 


+ W{0,m)]. (45) 


0) = ^([(iVi,e + iVi,o) - {N2,e + N2,o)f)- (48) 

Since hopping is permitted only along the chain, the 
particles will distribute themselves equally between the 
chains in order to minimize the energy. Hence, for the 
full span of g 


5 '(7r,0) = 0, 


(49) 
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a fact verified by the depicted values of S{Tr,0) in panel 
(c) of Figs. EHHiniini and[71 Additionally, we also have 
from Eq. (l45ll 

5(0,0)= ( ^ + ^ ] [lF(l,m) + lF(2,m)]. 

\77ieven ^odd/ 

Then, from Eq. (jC14ll in the Appendix [Cj we have 

mO) = ^([^p"-iVp']), (50) 

where Np = {Ni^e+^i,o) + {N 2 ,e+^ 2 ,o) is the total parti¬ 
cle number operator. Since we are working in a canonical 
ensemble and since the ground state is an eigenstate of 
Np, (at all couplings g) we get the simple relation 

5(0,0) = 0, (51) 


which is manifested in panel (d) of Figs. [31 SI O [51 and 

m 

To understand the phase transition and the ordered 
state, we will now study some aspects of the correlation 
function and the structure factor. We will consider an 
extreme situation of the ordered state, namely, only one 
of the sub-lattices of the two-chain system is occupied. 
Then, in Eq. (i43l) . for I = 1 (Z = 0) and m being even 
(odd), we have {nz2.j^k'nz^-,j+i,k+m) = 0- Subsequently, 
we obtain 


and 


47V2 

VF(0,modd) = -^, (52) 

4iV2 

VF(l,meve„) = -^. (53) 


When only one of the two sub-lattices is occupied, the 
ground state becomes an eigenstate of the operators 
Nj^o- Consequently, Eq. (H5|) and Eq. (ITT)) simply as fol¬ 
lows: 


4JV2 

[5(7r, 7r)]niax = (54) 

and 


[5(0,7r)]„,i„ = 0. (55) 

Based on the above considerations, we will analyze the 
general features in the correlation functions and structure 
factors at various filling fractions in the interacting CJT 
chains by studying the Figs. [31 [H El [H and [71 In all 
these figures, we consider the system at three coupling 
strengths: (a) when the system is quite disordered; (b) 
when the system is around the transition point; and (c) 
when the system is deep inside the ordered phase. 

First, panel (a) in all these figures shows the two-point 
inter-chain correlation function As the inter¬ 

action strength is increased, the curves become more os¬ 
cillatory; at large coupling, W{1, m) attains the constant 


4Af 

minimum value — at all even values of m as predicted 
by Eq. (1531) for the completely ordered case. Second, 
panel (b) in all the Figs. [31-[71 portrays the variation of the 
two-point intra-chain correlation function W (0, m). Here 
too, the oscillations become more pronounced as the cou¬ 
pling strength g increases. Furthermore, at very strong 
coupling, W{Q,m) attains the same minimum value of 

[as given by Eq. (1531) ]. but at odd values of m in 
contrast to W{l,m). 

Next, we observe that there are only two allowed val¬ 
ues for the wavevector component ky, i.e, it and 0. We 
study the behavior of the structure factors S{'K,kz) and 
5(0, kz) in the panels (c) and (d) of Figs. [5]-[7| In panel 
(c) of these figures, at wave vector (7r,7r), 5(7r, fc^) grad- 

4Af^ 

ually attains its maximum value of —[in agreement 
with Eq. m] as coupling becomes stronger. It is in¬ 
teresting to note that, in the disordered phase (i.e., at 
small coupling) 5(7r, k^) has a peak at kz = ^ttv which is 
reminiscent of the peaks in the Luttinger liquid for the 
CBM model (of Ref. and for the Holstein model (of 
Ref. i); furthermore, the peak at k^ = 2tiv becomes less 
pronounced as the coupling g increases in our interacting 
chains. This peak corresponds to only short-range corre¬ 
lations; consequently, it is expected to be not observable 
in a thermodynamic system. Last, in panel (d) of all the 
above-mentioned figures, 5(0, tt) gradually decreases as 
the interaction strength increases and finally drops to its 
lowest value of 0 as estimated by Eq. (1531) . 

We will now define the order parameter 5* (tt, tt) as the 
following rescaled value of 5(7r,7r): 


a*i 5(7r,7r) - [5(7r,7r)], 

5 (7r,7r) = 


[5(7r,7r)]niax “ [5(7r,7r)]i: 


(56) 


where [5(7r, 7r)]inin and [5(7r, 7r)]inax are, respectively, the 
minimum and the maximum values of 5(7r,7r). Hence, 
upon tuning electron-phonon coupling g at any general 
filling up to half-filling, the order parameter 5*(7r,7r) 
varies from zero to one when the system transits from 
a charge disordered phase to a CDW state [with order¬ 
ing wavevector (tt, tt)] as displayed in Fig [5] Interest¬ 
ingly, similar to the CBM model of Ref. [^, our interact¬ 
ing CJT model also predicts a conducting commensurate 
CDW state without an excitation gap . 


B. Ground-state fidelity and fidelity susceptibility 

Although we analyzed the phase transition using the 
rescaled structure factor 5*(7r,7r) in the previous sub¬ 
section, to pinpoint the transition better we will use the 
ground-state fidelity and the fidelity susceptibility. The 
ground-state fidelity, i.e., the overlap of the ground-state 
wavefunction at two very close values of the control pa¬ 
rameter (i.e., at g and g + S), is defined as follows^: 

F{g,6) = |(4'o(5)|4'o(5 + <5))|, 


(57) 
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g 

FIG. 9. (Color online) Plots of (a) ground-state fidelity and 
(b) fidelity susceptibility evaluated at <5 = 0.05 and at adia- 
baticity t/cjo — 0.1. 


eracy after phase transition, we invoke the same method 
devised in our earlier worbi^. 

In Figs. Ha) and Hb), we display our calculated 
ground-state fidelity and the corresponding fidelity sus¬ 
ceptibility for interacting-CJT-chain systems of different 
sizes and at various values of strong-coupling. The dips 
in the fidelity and the concomitant peaks in the suscep¬ 
tibility are indicative of the location of the phase transi¬ 
tion in the ^-parameter space; these transition points are 
in agreement with the order parameter curves in Fig. [8] 
Hence, ground-state fidelity and fidelity susceptibility are 
useful tools for understanding phase transition in our in¬ 
teracting CJT model. 

In the following section, we shall study interacting CJT 
chains in the adiabatic regime and at weak coupling. 

V. WEAK COUPLING ANALYSIS OF 
INTERACTING IDENTICALLY-LONG CJT 
CHAINS 

To analyze the weak-coupling case, we introduce new 
notation for ease of manipulation. We adopt the no¬ 
tation that, for the (ij ,2 orbital located at site (j, k) in 
Fig. HD c) and denoted by the destruction operator 
the surrounding oxygens in the y-direction are located at 
Uy.j_i j^ and while the adjacent oxygens in the 

z-direction are located at and Conse¬ 

quently, the corresponding second-quantized operators, 
analogous to those in Eq. ([34]), are given as 


where jlko) is the ground state of the system and S is 
a small increment in g. Next, fidelity susceptibility^ 
is defined as the second derivative of the ground-state 
fidelity: 

XF(ff) = = 21ini (58) 

< 5->.0 0 

It should be pointed out that, while the fidelity is depen¬ 
dent on the chosen value of 5, xpig) is independent of 
it. For numerically evaluating the values of fidelity and 
fidelity susceptibility, to avoid the problem due to degen- 

I 


Qx‘,j,k 


'\/2A(fc<Jo 


2fwo 


(59) 


Then, for the interacting CJT chains, we rewrite Eq. dSSj) 
as follows: 



+ H-C.) + Wo 


i.fc 

3,k 


j,k 

+az-j,k+^) - + ^z-,j,k-k 

)} + 4{(^y;j+l,fc + by-j+l,k) “ 

1 /1 

+ 2 v/2 i^z^-j,k- 

(60) 


Next, we Fourier transform the above Hamiltonian using 
the following: 

d 2 i. — ^ J(.P3+qk)g 

p,q 




— ^ pi[p3+q{k+^)] 

p,q 


'^\p,qi 
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and 




1 

7 ^ 


/ , ® ^y,p,q^ 

p,q 


r ■ i. — ^ \^^i(pi+qk) 

^x-],k — /Jr ^a:;p,5- 

p,q 

We express where 

ttWC _ _ \ J 

^0 — / . ^P,q^z^:p.q^^^;P,Q 

P,Q 

^ ^y;p,q^y,P,Q ^x;p,q^x-,p,q) i 

P,Q 

(61) 

with Cp^q = —2t cos(g) and 

^)Pp^q^ (62) 

^ P,<3 

with Pp^q = J2p>,q/ dl 2 .p,+p^q,+qd, 2 .p,^q, and 

QeffiP^ P) = 2* sin ^ (a^._p + a^-p^q) 

+ ^sin(|) (4;-p.-g + ^y;p.9) 

+ ^'^^(7;-p,-(j + '^aiip.q)- (63) 

We note that Q^f / (p, g) is proportional to the Fourier 
transform of v^Qi + Qa where Qi and Qs are the breath¬ 
ing mode and a Jahn-Teller mode^, respectively, at site 
(j, fe). The double time derivative of Qeff{p,q) is given 
by 

Qeff{p,q) = - [[Qeff{p,q),y^Y . (64) 

Let the eigenstates and eigenenergies of be denoted 
by \ipi) = \n)ei 0 \mx,my,m^)ph and respectively, 

with |0o) = |0; Oa;, Oy, Oz) being the ground state with 

zero phonons. Furthermore, and are the corre¬ 
sponding eigenstates and eigenenergies, respectively, of 
the interacting Hamiltonian + El^^ ■ 


Then, the matrix elements of the above Eq. (l64l) can be 
written as 

{^i\Qeff{p,q)\‘^i') = - {E^i - E/^pY ($i|(5e//(p, (?)|$p). 

(65) 

From the above equation, we see that during the tran¬ 
sition from |4>//) to |4>/), instability occurs if = 

- E^pY < 0 provided that {^i\Qeff{p,q)\^i') ^ 0. 
In the weak-coupling regime, we have 

^^7/ = ~ = {^o + g'^i^o Rexo(p,g;wo)}^ 

Ri Wq {l -f 27^0 R.exo(p,9;wo)} , 
where Xo{p^ 9! wq) is the relevant non-interacting dynamic 
polarizability defined below in Eq. (ED. Therefore, the 
instability condition is given by 

1+2g^uJoRexo(,P,q;^o) =0- (66) 


For a detailed and rigorous derivation of this instability 
condition (involving the real part of the non-interacting 
dynamic susceptibility), we refer the reader to Ref. @. 

Now, up to second order in perturbation, the energy 
of the ground state of is given by 


E,s>o = 


= Ko- 


ICJT 


- Y 

1^0 


E° - El 

4>i <Po 


2 2 
g ^0 

N 


Y 


p,q,m^0 
\{m\elPp,q\0)el\'^ 
^mO + 


/ 4 sin^ (; 

[) -f i sin^ 1 

f?) 

' + 5 } 

1 

iJ 4 

\2J 

8J 


(67) 


where ^mo = — Co with Cm being the energy of the 

state \m)ei. 

Next, we consider an excited state of EIq'^ with a single 
phonon of momentum (—p, —q) given by the general state 


|'0°) = /3a|0; Oa:, 0^, lz;_p,_,) -I- /3\/l - a^jO; O^,, ly-p-q, 0^) 
+Vl-/32|0;l.;-p.-g,0j„0z), (68) 


where — 1 < a < 1 and — 1 < /3 < 1. We evaluate the 

eigenenergy E^o for the interacting state 14'°) of 

as: 


E^o = {fi\B m)- 2 ^ - ^ - 




- E^q + Wo - 


2 2 
g ^0 

N 


Y. {'‘™'(|) + ^i”'(|) + 

r/ ^-An V 


1 1 |(m|eiPp,g|0)ei| 


p,q,m^0 


8 J CmO + Wo 


N 


Y 

m^O 


2/3a sin -f i/37l - ^ sin ^ 


I (^|e//^p,(j |0)e/1 I ^7Tz| 10) e/1 


^mO H" ^0 


CmO ^0 


( 69 ) 
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After subtracting Eq. (l67l) from Eq. (|6^ . we obtain 

= Wo + Re xo(-p, -<?; wo), 

where 


(70) 


Xo(p,g;wo) = |(^2/3asin (I) + i/3\/r^^7r2sin ^ + ^(1 -/3^)|xo(7',g; wq), (71) 


with 


1 


Xo(p,g;wo) = 


N 

m^O 

Based on the derivation in Appendix [D1 we have 

1 


\{m\elPp,q\0)el\'^ _ \{m\elP-p-q\0)el\ 

Wo - Cmo +ir] Wo + CmO + ip 


(72) 


Rexo(7',g;wo) = 


7 


27rwo Y^l - 72 


In 


(l + Vl-72) - {qtan(%^ + f)} 
(l-0^)^-{7tan(^ + f)}^ 


— In 


(l + yW)^-{7tan(-^ + |)}- 
(l-yW)^-{7tan(-^ + |)} = 


(73) 


where 7 = < 1- For values of the adiabaticity parameter 2 < t/ojQ < 5 that occur in manganites^ and for the 

range of filling fractions 0.05 < v < 0.5 that include C-type antiferromagnetism in manganites, it should be pointed 
out that 7 < 1. In the above expression, it is interesting to note that Xo(Pi9;wo) is independent of p. To obtain the 
critical coupling gc at which the instability first sets in, we need to maximize Xo(Pj9;wo) with respect to p. It can 
be shown that the coefficient {2/3asin (|) + |/?\/l — sin (|)} + |(1 — (3^) [occurring in Eq. (ED] is maximized 

4 sin(-) 

for p = TT with a = , and /3 = 1. This value of p = tt is expected because of inter-chain particle repulsion. 

\/l + 16sin2(f ) 

Therefore, the maximum Xo(7')'?jWo) is given by 


[Rexo(7r,g;wo)]j- 


= -—5-|i 

27r 4wo ( 


-1-16 sin' 


— In 


(!) 




In 


(l-bVl-72) - { 7 tan(%^-f f)} 
(1-Vl-7^)^- {7tan(^-f f)}^_ 


(l + yr^)^-{7tan(-^ + f)}- 
(l_^/r3^)2_{7tan(-^ + f)}^ 


(74) 


Using the instability condition 1 -|- 2p^wo [Re Xo('^i 9! Wo)]max = 0 with q = 2kp^, we get the critical coupling for 


filling fraction v = kp / ‘<is 

;{ 


^ = i|l-I-16sin^(A:i?)| 


5, 


^1-7^ 


In 


(1 - - 72) - (7tan(fcF))" 

(1 -b y/l - 7^)^ - (7tan(fcF))" 


- 21n 


1 - 1 - 7^ 

1 + Vl - 7^ 


(75) 


Based on the above Eq. ED, Fig. [ini displays the small 
parameter values at the transition point g = gc for var¬ 
ious fillings (up to half filling) and for different values 
of the adiabaticity parameter. Our perturbation theory 
is valid when the small parameter pcWo/t < 1- For the 
critical couplings pc in Fig- UHl the system undergoes a 
quantum phase transition from a disordered state to an 
insulating CDW state with ordering wave vector (tt, 

As t/uJo increases, even for smaller filling fractions, phase 
transitions are predicted within our approach. Interest¬ 
ingly, for values of the adiabaticity parameter 2 < t/ujQ < 
5 that are relevant to manganites^. Fig. [Tni portrays the 


critical coupling for density-dependent CDW ordering at 
filling fractions where C-type antiferromagnetism is real¬ 
ized in manganites^^— . 


VI. DISCUSSION 

We will now present, at strong coupling, a simplifi¬ 
cation of the model for the interacting two-chain CJT 
system [represented by Eq. (H^ j by mapping it onto a 
model for a single CJT chain. For ease of discussion, for 
the two-chain case, we first recast Eq. with general 
















































15 



V 


FIG. 10. (Color online) Plot of small parameter at critical 
coupling (leading to density-dependent CDW transition) for 
various filling fractions v and adiabaticities t/ojo. 


coefficients: 

-^e// =~^l + H.c.) 

j,k 

J._]^(l — 2n2;2y fc) 

X (1 + ffi H.C.] 

+nE nz‘^-j,knz^-j,k+i 
o,k 

+V'2^n^2.j^fen22y+i_fc, (76) 

j,k 

where Vi —>■ oo and V 2 —>• oo. The two chains can be 
superposed on each other (to form a single new chain) 
with the following conditions: (a) inter-chain particle re¬ 
pulsion V 2 —>■ 00 implies that in the new chain there is 
also hard-core repulsion between the particles originating 
from different chains; (b) to have more number of un¬ 
blocked sites for particles to hop to (so as to minimize the 
energy), particles originating from different chains should 
alternate along the propagating z-axis of the new chain; 
and (c) because Vi —>■ 00 , the next-nearest-neighbor hop¬ 
ping for particles originating from both the chains has the 
modified form -t 2 d^ 2 .fe_,_i(l-n^ 2 .fc)(i^ 2 .fe_i where is 
the creation operator of a particle at site k of the new 
chain and riz^-k = On incorporating the above 

conditions, the of Eq. (1751) has a simplified mapping 
onto the following Hamiltonian of a single chain: 

Hf.ff = —ti + H.c.) 

k 

“^2 -b H.C.]. (77) 

k 


Thus a system of N sites and Np particles in the 
interacting-chain model (given by has been re¬ 

duced to a system of iV/2 sites with Np particles [and 
governed by Eq. GZD], thereby leading to a significant in¬ 
crease in the size of computationally accessible systems; 
the number of basis states in the explored Hilbert space 
reduces from ^Cnp to "^CNp. Furthermore, the Hamil¬ 
tonian for a single CJT chain as given by Eq. (EH), for 
a system of N sites and Np particles, can be reexpressed 
in a generalized form 


^eff = -h + H.c.) 

k 

-t2 ^ ~ ^nz'^-,k)dz^-k+l + H.C. 

k 

+y'^'n^2.knz'i-k+i, (78) 

k 

with H > 00 . Then, as pointed out in Ref. M the model 
of the above Eq. ((751) can be mapped onto the following 
model with N — Np sites and Np particles: 


H. 


RC 

eff 


-—h ^(dl 2 .^d^ 2 .,, + H.c.) 

k 

-t 2 ^ “ n^ 2 .f,)d^ 2 .^ + H.c. 


(79) 


Here, it should be noted that the Hamiltonians in Eqs. 
dza) and dZH are identical. From the above considera¬ 
tions, it follows that the evolution of the order parameter, 
fidelity, and fidelity susceptibility in the interacting two- 
chain system with N = n sites and Np = Up particles 
[and represented by Eq. (1171) ] is identically mimicked 
by a single chain system with N = Up + n/2 sites and 
Np = Up particles [and represented by Eq. (1751) ] as shown 
in Figs. 1TT1 and 1T71 for n = 16 and Up = 4. Additionally, 
we also note that the nature of the phase transition for 
the CBM model of Eq. (1751) [and studied in Ref. [l^ 
is the same as that for the one-dimensional CJT system 
represented by Eq. (1751) because of the close similarity 
of the two governing equations and the structure factors 
depicted in Fig. [7J Furthermore, since the CBM model 
of Ref. [H undergoes a second-order phase transition and 
because the interacting-chain JT system can be mapped 
onto the single-chain JT system, we expect the phase 
transition (that occurs when the coupling g is varied) in 
interacting CJT chains to be also second-order. 

Next, we would like to present the relevance of our 
derived model (for interacting CJT chains in the strong 
coupling regime) to a different system. A wide variety of 
condensed matter models can be simulated by suitable 
engineering of ultra-cold atom systems in optical lattices. 
We would like to point out that our studied model of in¬ 
teracting chains can be realized in a cold-atom system 
represented by hard-core-bosons (HCBs). To this end, 
we transform a fermion (represented by a creation oper¬ 
ator d^a) to a HCB (represented by a creation operator 
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FIG. 11. (Color online) Plots of the order parameters veri¬ 
fying that a system of interacting CJT chains with = 16 
sites and Np = 4 particles is closely mimicked by a single CJT 
chain with N = 12 sites and Np = 4 particles. 



g 

FIG. 12. (Color online) Plots of (a) ground-state fidelity and 
(b) fidelity susceptibility showing very similar phase transi¬ 
tion for a system of interacting CJT chains (with = 16 
sites and Np = 4 particles) and a system of a single CJT 
chain (with A'' = 12 sites and Np = 4 particles). 


h'^) by using the relation h] - - 2 n^ 2 .j j,); 

this transformation is equivalent to a Wigner-Jordan 
transformation followed by a spin-to-HCB transforma¬ 
tion. Then the equivalent of the fermionic model in Eq. 


(1421) is given by the following model for HCBs: 

tt2D 

^HCB 

= —ti 'y -|- H.c.) 

-i2 - ii^+i,k)hj,k+i + H.C.] 

j,k 

+^1 ^Ik^lk+I + ^2 X! ^Ik-rij+i^k, (80) 

j,k j,k 

where (= hj Jij^k) is the number operator for HCBs. 

As pointed out in Refs. andH^, nearest-neighbor and 
next-nearest-neighbor hopping terms can be tuned inde¬ 
pendently; additionally, nearest-neighbor repulsion can 
also be achieved^. Thus, ladders and two-dimensional 
systems in general [which realize the Hamiltonian in Eq. 
(1801) ] can be simulated experimentally for a broad range 
of values of the parameters fi, ^ 2 , Vi, and V 2 • It is in¬ 
teresting to note that a one-dimensional version of the 
above Hamiltonian, i.e, 

^HCB =-ti “ ^2 J2(^k-i^k+i + H.c.) 

k k 

+ H^n^4+i. (81) 

k 

was studied by Mishra et a h^^d^ . These authors consid¬ 
ered kinetic frustration and obtained the phase diagram. 
Eurthermore, in the extreme case of ti = 0, the above 
model [of Eq. dSTll ] was recently shown to undergo a 
striking discontinuous transition from a superfluid to a 
supersolid in Ref. 


VII. CONCLUSIONS 

In this paper, we studied the CDW transition in in¬ 
teracting CJT chains at strong coupling [characterized 
byi^ t/{gujo) < 1] and at weak coupling [identified by^ 
{gujQ)/t < 1]. The nature of CDW in the interacting CJT 
chains is quite different in the two cases: strong coupling 
analysis predicts a conducting CDW with Z 2 symmetry 
being broken (i.e., one sub-lattice has higher density), 
whereas weak coupling theory shows that an insulating 
CDW results with a wavevector that varies linearly with 
density. Our study is relevant in manganite systems with 
C-type antiferromagnetism and suggests identifying the 
charge-ordering wavevector as an alternate way for deter¬ 
mining the controversial regime of the electron-phonon 
coupling. In the case where the coupling is strong, we 
use structure factor and fidelity to track the quantum 
phase transition. Exact expressions for the structure fac¬ 
tor at special wavevectors [i.e., (0,0), (7r,0), (0,7r), and 
(tt, tt)] are derived to understand the evolution of the sys¬ 
tem when the tuning parameter is varied. On the other 
hand, in the weak coupling regime, we identify the crit¬ 
ical coupling at which the phonon mode (involving the 
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breathing mode and a Jahn-Teller mode) becomes soft; 
our analysis is applicable for the adiabaticities occurring 
in manganites. Although our analysis for the interacting 
CJT chains corresponds to a two-dimensional case, it is 
extendable to three-dimensional systems as well and the 
conclusions are expected to be similar. 
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= - l)(e«' - l)(e9' - 1) + - l)(e«' - 1) 
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= - 1. (A3) 


We thank P. Majumdar, P. Littlewood, N. D. Mathur, 
D. E. Khmelnitskii, and A. Ghosh for useful discussions. 


Appendix A: Simplification of the function 

G7i(oi, 02, . . . , Oti) 


In this appendix, we obtain simple expressions for the 
function G„(ai, 02 ,..., a„) appearing in the main text. 
The general term G„(ai, 02 , ■■■,««) is defined as 


Gn(o^l,0^2; ■ • ■ jCTyi) — Ffi (oi, 02, . • . , Oin} 
n — 1 

+ X/ X/ ^kictci 1 ^C2 7 ■ * ■ ; ; 

k—1 c 


where 

00 

-^n ; ■ • ■ ; ^ ^ 

m\ — l 


E 


rrin — l 


(oig^)"^! ■ ■ ■ (ang^)”"" 
mi!... m„!(mi m„) ’ 


and the summation over c represents summing over all 
possible "Gm combinations of m arguments chosen from 
the total set of n arguments {oi, 02 ,..., an}- 

We begin by examining the simple case of the term 
G2(2,2) [= E2(2,2)-k2Fi(2)] appearing in Eq. We 

evaluate the derivative, with respect to g^, of G 2 ( 2 , 2) as 
follows: 

g^^G2i2,2)= - 1) + - 1) 

= {(e2s'-l) + l}{(e29'-l) + l}-l 
= - 1. (Al) 


Then, on performing integration, we get 

- 1 


G2(2,2) = J ■ 




-dg 


- 


n=l 

^ ^2\n 


= E 




n n\ 


(A2) 


Next, we consider the other term G3(4,1,1) [= 

^3(4,1,1) + 2 ^ 2 ( 4 ,1) + ^2(1,1) + El(4) + 2Ei(l)] occur¬ 
ring in Eq. (EZl). On taking the derivative of G3(4,1,1) 


Hence, we get 

G3(4,1,1) = J ■ 


, 69 = _ 1 


-dg 




n—1 


n\ 


= E 


( 65 ^)^ 


d ■ 


(A4) 


Finally, we evaluate the derivative of the general term 
G„( 01 , 02 , ■ • ■ ,Q:„) with respect to g^: 

g ^ ^G7i(oi , O 2 , . . . , On) 

= (e“is' - l)(e“"«' - 1)... (e“'‘9" - 1) 

n—1 

+ ^ - l)(e““=9' - 1)... - 1) 

- 1 


k—1 c 


= [n(Li{(e“-'^'-i) + i} 

= qT,7=1 _ 1, 

Then, the general term is obtained to be 


(A5) 


Gn(oi, 02, . . . , 0„) — 


gELi “is" _ 1 






m—l 
00 vn 


= E 


(Er=i«.<?^) 


ml 

2\m 


m ml 


(A6) 


For large values of g^, we have the approximation 


5Er=i °‘i9^ — 1 


-dg^ 


eEr=i “iS 

E n 0 ' 

i=i 


Appendix B: Comparison of phase transition in 
interacting CJT chains at strong coupling and 
various values of adiabaticity 

In this appendix, we compare the behavior of the in¬ 
teracting CJT chains at strong coupling and for different 
values of the adiabaticity parameter t/uJo- In Fig. [13] 
we portray the order parameter S'*(7r,7r) for two fill¬ 
ing fractions and for two sufficiently different values of 
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FIG. 13. (Color online) Plots comparing order parameter 
S*{'K, tt) at different values of t/ujo and at different filling frac¬ 
tions. 


Appendix C: Exact expressions for the structure 
factor in limiting cases 

In this appendix, we derive formulae for S{'!r,Tr), 
5(0, tt), 5(7r, 0), and 5(0,0). The structure factor is ex¬ 
pressed in terms of the correlation function W{l^m) as 
follows: 

S{ky, k,) = m), (Cl) 

l,m 


where ky = Uyir with Uy = 1,2 and kz = with Uz = 
1,2,...,-^. Since we are dealing with two interacting 
chains, the above equation simplifies to 



-hlE(0,m)].(C2) 


From the above Eq. (IC2I) . it follows that the structure 
factor at (tt, tt) is given by 



g 

FIG. 14. (Color online) Plots comparing (a) ground-state 
fidelity and (b) fidelity susceptibility at 5 = 0.05 and Np = 4 
when using different values of t/ujo and fillings Np/N. 


N/2 


5(^, tt) = ^ (-1)™[-1E(1, m) + 1F(0, m)] 

m—1 

= ( E - E ) [-^(l’H + m0,m)], (C3) 

\77leven ^odd/ 


at (0, tt) is given by 


N/2 


5(0, tt) = ^ (-1)-[1E(1, m) + W{0, m)] 

m—1 

= E-E [W{l,m) + W{0,m)], (C4) 

\?71even ^odd/ 


at (tt, 0 ) is given by 


N/2 


5 (7r,0) = ^ [—FF(l,m) -|- FF(0,m)] 

m—1 

= ( E + E) + (C5) 

\?7leven ^odd/ 

and at (0,0) is given by 


the adiabaticity parameter. As is evident from the fig¬ 
ure, for both the fillings, the transition for t/oJo = 0.5 
takes place at smaller values of g compared to the transi¬ 
tions for t/uJo = 0.1. This observation is consistent with 
the transitions indicated by the extrema in the plots of 
ground-state fidelity and fidelity susceptibility displayed 
in Fig. [TH 


N/2 

5(0, 0) = ^ [W{1, m) + 1F(0, m)] 

m—1 

E+E [VF(l,TO)-bVF(0,TO)]. (C6) 

^even ^odd/ 

Now, we calculate the various terms in the above equa- 
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tions as follows: 


Next, we observe that 


E 


m) 


^ W{Q,m) 

^even 


= -y 

N 


N‘^ 

/ \ p 

;j ,k'^ 


even j,k L 

’T-z^y.fe+m) ~ E/ 7^ 


j,k 




= -E 

/V 


^ 2 ;^ ;j,fc+m) 


2^ ;_7,fc+m) EE 




= -E 

TV ^ 


= - E E 

TV ^ ^ 


^even 


TV^ 

\^z^]j,k‘^z^-J-\-l,k-\-m) j^2 


= -E 

/V 2-^ 


T,fe 


= -E 

/V 


{klz^-j^k ^ ^ '^z2;j + l,fc+m) ^ ^ 

^even ^ever 


TV^ 

-i y p 

JP 


^2:^;j;-|-l,A:+m) EE 

Alodd ?^even k ^TZ-even 


TV^ 

i. y p 

Ip 


= -E 

TV ^ 


^ ^ (^z2;j,fc-^j + l,e) 
k 

'^even 

E ^ ^ ('^2^;j;,fc-^j + l,o) 
^odd 


±y p 


TTZ-even 

fe ?71even 

Tvf 

= -E 

TV 21 ^ 

I i 

’ . . . . TV^' 

(TV,.eTV, + i,e) + (TV,.oTV, + i,„) --f 

E 


_ p_ 

8 

0 

-fceven 

^odd 


4 



N 


2(TVi,eTV2.e) + 2 (TVi,„TV2.o) - -f 


(C8) 


Similarly, we also obtain 


^ bb(0,m) 


= -E 


N 

4 

TV 


TV^ 

(TV,,eTV,,e) + (AT,.oTV,.„) - ^ 


TV2 

2 \ P 


(Ee) + (Eo) + (Ee) + (Eo)- 


E ^(0’ 


where TV,-e = (^1.° = Efc^^d is the 

number operator for the total number of particles occur¬ 
ring in the even (odd) sublattice of the j-th chain. Then, 
it follows that 


= - E E 


’^odd j,k 


TV^ 

Alp 


{tI z^ ;j^k'^z^ -J,k-\-m} ^2 


= -E 

N Z^ 


j,k 


= -E 

TV ^ 


{Plz^-J^k ^ ^ '^z'^;j,k-\-m) ^ ^ 

, ^^odd ^^odd 

E E \j^k+m ) 


TV2 

P 

TV2 


TVi 


+Ei'^'^^b.fc EE E A72 


= -E 

TV 2 -^ 

j 

= -E 

/V 


(C7) 


TV 

4 

TV 


TV2 

/Codfj TTiodd fc TTlodd 

~ TVp 

/ . (^z2;,-,feTV,-,o) + / , {fkz^;j,k^j,e) ^ 

■^even ^odd 

' . . . . TV^ 

(TV,,eTV,.„) + (TV,.„TV,.e) - ^ 

TV2 


2(TVi,,TVi,„) + 2(TV2.eTV2,„) - 


(C9) 
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Lastly, we get 

W{l,m) 


= ^ 1 : 1 : 


^odd j',k 




= ly 

/V ^ 








ly 

N ^ 

E 

(n,22 

k ^ ^ '^z'^\j+l,k+m) 

3 

-^even 


^^odd 

+ y] 

E 


^odd 


^^odd 

k 

ly 

N ^ 

3 

E 

-feeven 

(71^2.,- 

fc^i+ 1 , 0 ) 


iV2 

i. Vp 

N2 


±y p 


= -y 

N 


N ■ 

4 

TV 


+ ^ ] {'^z^-j,k^j + l,e) Q 

kodd 

(TV,,,TV,+i,„) + (TV,,„TV,+i,,) - ^ 


TV 2 


2(TVi,eTV2,o) + 2(TVi,„TV2,e) - y 


(CIO) 


Using the expressions obtained in Eqs. (IC7l) . (ICSl) . S, 
and (jClOll . we simplify Eqs. (IC3I) . (IC4I) . (IC5E and 
as follows: 


<5(^,7r) = -([(TVi,, - TVi,„) + (TV2,„ - TV2,e)]'), (Cll) 


*5(0, tt) = _([(TVi,e-TVi,„) + (TV2,e-TV2,o)]2), (C12) 


*5(^,0) = l([(TVi,e +TVi,„) - (TV2,e + TV2 ,o)]2), (C13) 

and 

*5(0,0) = l([TVyTV2]), (CM) 

where Np = (TVi,e + TVi,o) + (TV 2 ,e + TV 2 ,o) is the number 
operator for the total number of particles. 
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Appendix D: Evaluation of dynamic susceptibility for interacting CJT chains at weak coupling 


In this appendix, we evaluate the real part of the dynamic susceptibility for interacting CJT chains in the weak- 
coupling regime. A simple analytic expression for the dynamic susceptibility is obtained as follows: 


1 


Rexo(p,<?;wo) = 


m^O 


\{m\elPp,q\0)el\'^ | (m|eiP_p 10)ei | 


=-y 

N ^ 


p .q L 


1^0 ~ ^mO Wo + ^mO 

(1 +p^q'^q' ^z'^;p' —p,q' — q^'^z'^;p\q' 


=-y 

AT 


N 

1 


p',q' L 

f'kp 


•^0 ^p',q') 

*^0 “ 1 “ i^p' —p,q' —q ^p',q 

) _ 

’kiz^;p' ,q' 

3 


Wo ~ {^p'+p,q'+q ~ ^p',q') 

^0 {^p'-\-p,q'-\-q ^p',q') 


-1 

1 

1 



-kp 

pkp 


Wo — 4t sin(q' + |) sin(|) Wq + 4t sin(g' + |) sin(|) 


J Wo 


dq' 


dq' 


7 -sin(g' + |)) 7 + sin(g' + |)) 


(Dl) 


where 7 = 4 ^ g ^ and use has been made of the facts that there is reflection symmetry and that ep^q is independent 
of the value of p because there is no inter-chain transport. Integrating the above equation for the case 7 < 1, we 
obtain 


R.eXo(p,9;wo) = 


1 


7 


27rWo y^l — 72 
-In 


In 


(l + vT^)^-{7tan(^ + f)}^ 
(l-vT^)"-{7tan(^ + f)}^ 
(l + yr^)^-{7 tan(-^ + |)}^' 


(1-yr^) -{7tan(-^ + |)}^ 

I- 


(D2) 
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Appendix E: Popular Summary 

Understanding the exotic phenomena (such as colossal 
magentoresistance, coexisting charge-, spin-, and orbital- 
orderings) in bulk manganites and designing artificial 
structures (such as heterostructures, quantum wires, and 
quantum dots) using manganites is of immense funda¬ 
mental interest and also of huge technological importance 
(especially for development of electronic and spintronic 
devices). 

To study emergent ordering, exploit device potential, 
and guide material synthesis in these complex magnetic 
oxides, effective Hamiltonians are needed for various in¬ 
teractions. Except for the cooperative electron-phonon 
interaction (EPI), effective Hamiltonians (that reason¬ 
ably model the essential physics) have been derived for 
all other interactions. For instance, large Hunds cou¬ 
pling can be mimicked by the double exchange model, 
interactions between localized spins (that have strong 
on-site repulsion) can be described by superexchange 
model, and Hubbard Coulombic interaction can be mod¬ 
elled by dynamical mean-filed theory. The main chal¬ 
lenges for modelling cooperative EPI in manganites have 
been mathematical analysis of the quantum-phonon ef¬ 
fects and resolving the ambiguity in the nature of the 
electron-phonon interaction. Here, using a controlled an¬ 
alytic treatment of quantum phonons, we show that weak 
coupling and strong coupling produce different experi¬ 
mentally observable effects. 

Since manganites with C-type antiferromaganetic or¬ 
dering can be effectively considered as chains, we study 
the effect of electron-phonon coupling in interacting 
chains. At strong coupling, our theory predicts a con¬ 
ducting charge-density-wave (CDW) with charge resid¬ 
ing predominantly in one sub-lattice; hence, the order¬ 
ing wavevector is independent of density. On the other 
hand, at weak coupling, our theory shows that an insulat¬ 
ing CDW results with a wavevector that varies linearly 
with density. Thus, based on the experimental deter¬ 
mination of the charge-ordering wavevectors (at various 
doping values) in manganites that exhibit C-type anti¬ 
ferromagnetism, we believe that our theory offers an op¬ 
portunity for identifying the regime of electron-phonon 
coupling. 


